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PRELIMINARY AMENDMENT 

Assistant Commissioner of Patents December 1 8, 2001 

Washington, D.C. 20231 

Sir: 

Prior to examination of the above-identified application, please amend the 
application as follows: 

IN THE SPECIFICATION 

Please replace the original specification with the following substitute 
specification, as attached hereto. A marked-up version of the original specification 
identifying the amendments in the Substitute Specification is also attached. 

IN THE CLAIMS 

Please amend claims 1-7, and add new claims 8-13 in this new application. A 
marked-up version of the amended claims will follow on proceeding pages. 
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(Amended) A method for determining reflection travel times of seismic 
events picked on 3D records of seismic data corresponding to signals 
picked up by seismic receivers distributed along an acquisition line in 
response to the emission in the subsoil of waves from source points, the 
seismic data having been previously converted to cylindrical wave data 
and interpreted, comprising: 

a) defining a slowness vector (p ) whose component (p x ) in a direction 
parallel to the acquisition line defines the slope of the cylindrical 
wave; 

b) for a seismic receiver situated in an abscissa (X R ) on the acquisition 
line, seeking an abscissa (f) of the source point on the acquisition 
line such that a ray starting from the seismic receiver and reflecting 
on a picked event emerges at the source point, with a slowness 
vector (p ) whose component in the direction of the acquisition line 

is (-p x ); 

c) determining a travel time (t e (X R )) by adding to the value of the travel 
time along the ray a time equal to the product of the slope of the 
cylindrical wave by the abscissa of the source point; 

d) repeating steps b) and c) for all the position of the seismic receivers 
for which a demigration result is wanted; and 

e) repeating steps a) to d) for all the values taken by parameter (p x ) 
for which an event has been picked in the corresponding migrated 
data. 



(Amended) A method for determining reflection travel times of seismic 
events picked on 3D records of seismic data corresponding to signals 
picked up by seismic receivers distributed along several acquisition lines 
in response to the emission in the subsoil of waves from one or more 
source points, the seismic data having been previously converted to a 
superposition of migrated cylindrical wave data and interpreted, 
comprising: 

a) defining a slowness vector (p ) whose component (p x ) in a direction 
parallel to the acquisition lines defines the slope of the cylindrical 
waves associated with the various acquisition lines; 

b) for a given seismic receiver situated at an abscissa (X R ) on an 
acquisition line, seeking an abscissa (f) of a source point on the 
acquisition line such that a ray starting from the seismic receiver 
and reflecting on a picked event emerges at the source point, with a 
slowness vector (p ) whose component in the direction of the 
acquisition line is (p x ); 

c) determining a travel time (t e (X R )) by adding to the value of the travel 
time along the ray a time equal to the product of the slope of the 
cylindrical wave by the abscissa of the source point; 

d) repeating steps (a to d) for all the acquisition lines for which a 
demigration result is wanted; 

e) repeating steps (a to d) for all the acquisition lines for which a 
demigration result is wanted; and 

f) repeating steps (a to e) for all the values taken by parameter (p x ) 



for which an event has been picked in the corresponding migrated 
data. 

(Amended) A method as claimed in claim 1 , wherein the travel times 
associated with the cylindrical wave events are used to implement a 
cylindrical-wave reflection tomography technique and to calculate a 
velocity distribution in the environment. 

(Amended) A method as claimed in claim 1 , wherein the travel times 
associated with the cylindrical wave events are converted to travel 
times associated with shotpoint events, by exploiting correspondences 
between the cylindrical wave data and the shotpoint data. 

(Amended) A method as claimed in claim 4, wherein the seismic data 
are organized by acquisition lines, the synthesized show travel times 
are converted to travel times associated with source-receiver pairs, the 
conversion being carried out by extrapolation, in the direction of the 
acquisition lines, of the travel time synthesized for a source-receiver 
pair that is closest to the source-receiver pair of the acquisition device 
for which information on the travel time is wanted. 

(Amended) A method as claimed in claim 4, wherein, the seismic 
data are constant-azimuth data, the synthesized shot travel times are 
converted to travel times associated with the source-receiver pairs, the 



conversion being carried out by means of a double extrapolation, the 
first one in the direction defined by the azimuth, the second in the 
direction orthogonal to the previous direction. 

7. (Amended) Application of the method as claimed in claim 4 to 
determine a velocity model by applying, a reflection tomography to 
travel times which are obtained. 

Please add new claims 8-13 as follows: 

8. A method as claimed in 2, wherein the travel times associated with the 
cylindrical wave events are used to implement a cylindrical-wave reflection 
tomography technique and to calculate a velocity distribution in the 
environment. 

9. A method as claimed in claim 2, wherein the travel times associated 
with the cylindrical wave events are converted to travel times associated 
with shotpoint events, by exploiting correspondences between the 
cylindrical wave data and the shotpoint data. 
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10. A method as claimed in claim 9, wherein the seismic data are organized 
by acquisition lines, the synthesized show travel times are converted to 
travel times associated with source-receiver pairs, the conversion being 
carried out by extrapolation, in the direction of the acquisition lines, of the 
travel time synthesized for a source-receiver pair that is closest to the 
source-receiver pair of the acquisition device for which information on the 
travel time is wanted. 

O 11. A method as claimed in claim 9, wherein, the seismic data are 

HI constant-azimuth data, the synthesized shot travel times are converted 

3f to travel times associated with the source-receiver pairs, the 

conversion being carried out by means of a double extrapolation, the 
j>3 first one in the direction defined by the azimuth, the second in the 

pi direction orthogonal to the previous direction. 



12. Application of the method as claimed in claim 5 to determine a velocity 
model by applying, a reflection tomography to travel times which are 
obtained. 



13. Application of the method as claimed in claim 6 to determine a velocity 
model by applying, a reflection tomography to travel times which are 
obtained. 
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REMARKS 

Claims 1-7 have been amended to remove the multiple dependent claims before 
filing fee calculation. New Claims 8-13 have been added. No new matter has been 
added to this application. 

To the extent necessary, Applicants petition for an extension of time under 37 
C.F.R. §1 .136. Please charge any shortage in fees due in connection with the filing of 
this paper, including extension of time fees, to Deposit Account No. 01-2135 (Case No. 
612.40916X00) and please credit any excess fees to such Deposit Account. 

Entry of the amendment is respectfully solicited. 

Respectfully submitted, 
ANTONELLI, TERRY, STOUT & KRAUS, LLP 




Donald E. Stout 
Registration No. 26,422 
(703)312-6600 

DES:clt 
Attachments 
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VERSION WITH MARKINGS SHOWING CHANGES MADE: 



1 (Amended) A method for determining reflection travel times of seismic 
events picked on 3D records of seismic data corresponding to signals picked 
up by seimic receivers distributed along an acquisition line in response to the 
emission in the subsoil of waves from source points, [this] the seismic data 
having been previously converted to cylindrical wave data and interpreted, 
[characterized in that it comprises the following stages] comprising : 

a) defining a slowness vector (p ) whose component (p x ) in a 

direction parallel to the acquisition line defines the slope of 
the cylindrical wave; 

b) for a seismic receiver situated in an abscissa (X R ) on the 
acquisition line, seeking an abscissa (f) of the source point 
on the acquisition line such that a ray starting from the 
seismic receiver and reflecting on a picked event emerges 
at the source point, with a slowness vector (p ) whose 
component in the direction of the acquisition line is (-p x ); 

c) determining a travel time (t e (X R )) by adding to the value of 
the travel time along [said] the ray a time equal to the 
product of the slope of the cylindrical wave by the abscissa 
of the source point; 

d) repeating [stages] steps b) and c) for all the position of the 
seimic receivers for which a demigration result is wanted; 
and 
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e) repeating [stages] steps a) to d) for all the values taken by 
parameter (p x ) for which an event has been picked in the 
corresponding migrated data. 

(Amended) A method for determining reflection travel times of seismic 
events picked on 3D records of seismic data corresponding to signals 
picked up by seismic receivers distributed along several acquisition 
lines in response to the emission in the subsoil of waves from one or 
more source points, [this] the seismic data having been previously 
converted to a superposition of migrated cylindrical wave data and 
interpreted, [characterized in that it comprises the following stages] 
comprising : 

a) defining a slowness vector (p ) whose component (p x ) in a 
direction parallel to the acquisition lines defines the slope of the 
cylindrical waves associated with the various acquisition lines; 

b) for a given seismic receiver situated at an abscissa (X R ) on an 
acquisition line, seeking an abscissa (f) of a source point on the 
acquisition line such that a ray starting from the seismic receiver 
and reflecting on a picked event emerges at [said] the source point, 
with a slowness vector (p ) whose component in the direction of the 
acquisition line is (p x ); 

c) determining a travel time (t e (X R )) by adding to the value of the travel 
time along [said] the ray a time equal to the product of the slope of 
the cylindrical wave by the abscissa of the source point; 



d) repeating [stages] steps a) to d) for all the acquisition lines for 
which a demigration result is wanted; 

e) repeating [stages] steps a) to d) for all the acquisition lines for 
which a demigration result is wantedM : and 

f) repeating [stages] steps a) to e) for all the values taken by 
parameter (p x ) for which an event has been picked in the 
corresponding migrated data. 

(Amended) A method as claimed in [any one of claims] claim 1 [or 2], 
[characterized in that] wherein the travel times associated with the 
cylindrical wave events are used to implement a cylindrical-wave 
reflection tomography technique and to calculate a velocity distribution 
in the environment. 

(Amended) A method as claimed in [any one of claims] claim 1 [or 2], 
[characterized in that] wherein the travel times associated with the 
cylindrical wave events are converted to travel times associated with 
shotpoint events, by exploiting correspondences between the 
cylindrical wave data and the shotpoint data. 

(Amended) A method as claimed in claim 4, [characterized in that] 
wherein the seismic data [being] are organized by acquisition lines, the 
synthesized show travel times are converted to travel times associated 
with source-receiver pairs, [this] the conversion being carried out by 
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extrapolation, in the direction of the acquisition lines, of the travel time 
synthesized for [the] a source-receiver pair that is [the] closest to the 
source-receiver pair of the acquisition device for which information on 
the travel time is wanted. 

(Amended) A method as claimed in claim 4, [characterized in that] 
wherein , the seismic data [being] are constant-azimuth data, the 
synthesized shot travel times are converted to travel times associated 
with the source-receiver pairs, [this] the conversion being carried out 
by means of a double extrapolation, the first one in the direction 
defined by the azimuth, the second in the direction orthogonal to the 
previous direction. 

(Amended) Application of the method as claimed in [any one of 
claims] claim 4 [to 6 for] to [determining the] determine a velocity 
model by applying, [for example], a reflection tomography to [the] travel 
times which are obtained. 
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VERSION SHOWING MARKED-UP SPECIFICATION 
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FIELD OF THE INVENTION 



The present invention relates to a method in^Scd to obtain Election travel times 
from an interpretation of seismic data in migrated cylindrical waves, for a given value 
of a parameter defining the slope of these waves, or from the superposition of such data 
5 associated with various substantially parallel acquisition line^this parameter possibly 
taking successively several valuesj 

BACKGROUND OF THE INVENTION-^ — ^ ^ — 

Demigration of seismic data is the invert of migration, which is a conventional 
seismic data imaging method essentially consisting, knowing the value of a wavefield at 
10 a known depth, for example at the surface, as well as a distribution model of the 
propagation velocities of the waves in the subsoil, in modelling the propagation of the 
source field and the backpropagation of the reflection data recorded, and in seeking 
phase coherences between these two modelled fields. Migration is particularly useful 
for interpretation of seismic data acquired on complex structures. 

15 It also allows to access, via demigration, the arrival times of the reflections 

associated with the picked events. This operation is notably carried out for determining, 
by means of kinematic methods such as migration velocity analysis methods, the 
velocity distribution in the subsurface. This determination constitutes a key stage in 
complex structures imaging. An example of such a method is described by : 

20 - Stork, C, 1992 ; « Reflection Tomography in the Postmigrated Domain » ; 
Geophysics, 57, 680-692. 
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A particularly effective migration velocity analysis method, referred to as Smart 
(« Sequential Migration Aided Reflection Tomography ») is notably described by : 

- Ehinger, A. and Lailly, P., 1995, Velocity Model Determination by the SMART 
Method, Part 1 : Theory : 65 th Ann. Internat Mtg., Soc. ExpL Geophys., Expanded 
5 Abstracts, pp. 739-742. 

Reflection tomography is used, in such methods, for determining the velocity 
mode!^;^e reflection travel times (reflection tomography data) are obtained, in the 
complex parts, by demigrating the picked events upon interpretation of the migrated 
data. This method has been successfully implemented for 3D imaging of complex 
10 structures by using 2D prestack migration algorithms (and the associated demigration) 
on a set of seismic profiles extracted from a 3D acquisition. However, such an approach 
is efficient only in a very specific context : 



on the geologic plane, the structtye must vary relatively slowly in a direction 



structure 

(referred to as « strike » mJEhe Angto^axoi^Uterature) ; 



15 - the seismic profiles must be acquired in the orthogonal direction (therefore the 
complexity direction) referred to as « dip » in the^n^o^^co^literature. 

On the other hand, this approach is difficult to use with complex structures without 
a privileged direction : such structures require prestack 3D migration techniques, which 
involve long calculation times, especially if imaging is carried out on a volume. 

20 \ ' Patent^R52,784,195^filed by the applieant^iotably describes another seismic data 
migration method knownfto the man skilled/ in the art as Alchemig. This method allows 
tef perform a 3D prestack migration of seismic events and imaging of the volumes of an 
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underground zone, with an attractive calculation time, from a series of a number N s of 
reflection survey cycles. 

It comprises successive emission of elementary wavefields, each defined by the 
association of a seismic signal W(t) and of a determined emission point from a series of 

5 emission points S : with 1 < i < reception, by seismic receivers placed at positions 

Rf , of the seismic signals reflected by the zone in response to each one of these 
wavefields, and recording of the various signals received by each seismic receiver in 
form of time-dependent seismic traces d}(f) , For a given velocity model, 

a) a slowness vector p (homogeneous to the inverse of a velocity) whose two 
10 components p y and p x can take a previously defined series of values is defined, 

b) a time lag function t 0 (pj) is defined for a given slowness vector p and for a given 
emission point S t , 

c) time lag function t 0 (pj) is applied to each elementary wavefield associated with 

emission point S { and a composite wavefield is formed at the surface by spatiotemporal 
15 superposition of the various elementary wavefields to which such a time lag is applied, 

d) a time lag t 0 (pj) is applied to each seismic trace df(t) marked by pair (i,j) and a 

composite trace field is formed at the surface by spatiotemporal superposition of the 
various seismic traces to which such a time lag is applied, 

e) a migration of the composite trace field is carried out by using as the wavefield the 
20 composite wavefield, which is done by modelling the propagation of the composite 

wavefield and the backpropagation of the composite trace field, and by suitably 
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combining the two composite fields thus modelled at any point of the zone to be 
imaged, 

f) sfcig^ c) to e) are repeated for all the values taken by components Py and p x of vector 
p , and 

5 g) for any fixed value of the second component p x of sarfi vector p , the results of these 
various combinations are summed so as to obtain a migrated image associated with this 
fixed value of p„ thus performing a prestack migration. 

It is assumed that, even if the technique proves robust in case of violation of these 
hypotheses, the seismic data is constant-azimuth data (the azimuth being defined as the 
10 direction of the source-pickup bipoints) (hypothesis 1) and that acquisition is carried out 
by displacing the source along lines parallel to the direction defined by the azimuth 
(hypothesis 2). 

If this is not the case, it can be remedied by first carrying out a preprocessing 
operation, for example according to the technique referred to as AMO (« Azimuth 
15 Moveout »), well-known to the man skilled in the art, to adjust the data obtained from 
standard marine acquisitions. 

The co-ordinates system selected here is such that the x-axis is parallel to the 
direction of the acquisition lines and the y-axis thus represents the position of an 
acquisition line^e defin^ Jtereafter^ wave datgX the^»^^^iflt, 

20 along an acquisition line, of the seismic response of the subsurface to a cylindrical wave 
(wave generated by a source line and whose phase varies linearly with the abscissa 
along the acquisition line) whose axis coincides with the acquisition line considered. 
Considering this definition, the Alchemig technique allows >6 obtain the superposition, 
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along the various acquisition lines, of /C migrated cylindrical wave data/, these 
cylindrical waves being associated with an acquisition line and with a predetermined 
value of parameter p x defining the slope of the cylindrical wave. The procedure consists 
in applying to the data linear phase shifts parametered by vector p whose component 
5 along the profiles axis is p*, and in calculating such a superposition by applying to the 
phase shifted data a plane-wave migration software. Postmigration stacking is obtained 
by summing the results obtained for the different values taken by p x . We suppose that, 
by suitable selection of the origin of the mark or of the phase shift function, the point of 
the surface where the phase shifts are zero has (x = 0, y = 0) as the co-ordinates. 

10 It can be noted that the technique described above is one way among many to obtain 

superposition of migrated cylindrical wave data associated with parallel acquisition 
lines and with a given value of parameter p x . Another procedure can for example 
-eerasist, as described by Claerbout, J.F., 1971, ^Towards a Unified Theory of Reflector 
Mapping,/ ; in Geophysics, 36, No.3, 467-481, in calculating the result by carrying out, 

15 for each acquisition line : 

- propagation of the cylindrical wave and backpropagation of the cylindrical wave 
data, and 

- calculation, at any point of the space, of the crosscorrelation of the source field 
whose propagation has been modelled and of the backpropagated data field, which 

20 allows J# obtain, for the acquisition line considered, the migrated cylindrical wave data ; 

and by stacking (superposing) the migrated cylindrical wave data obtained for the 
various acquisition lines. 
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If wcr-teve; for several values of parameter p*, the superposition of migrated 
cylindrical wave data^associated with various acquisition lines, it is possible, as 
explained in the aforementioned patent, to control the quality of the velocity model used 
for the migration and also to update the velocity model by exploiting the deformations 
5 of the events when going from one cylindrical wave parameter to the next (migration 
velocity analysis). 

For updating, according to the aforementioned Smart method (or others), the events 
bave to be picked, for various values of parameter p*, in the superposition of the 
Jn migrated cylindrical wave data , a superposition corresponding to the various acquisition 

m 10 lines, then the picked events have to be demigrated and the model is updated by 
reflection tomography. Even so one nofote to kno\^how to demigrate the interpreted 
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migrated events when interpretation has been performed on migrated cylindrical wave 
data or the|* superposition obtained for various acquisition lines. 

SUMMARY OF THE INVENTION 

15 £jhis is one of the complex problems that can be solved by the method according to 
the invention. The approach selected furthermore leads to methods that will be 
described, and which are closely related to the concept of cylindrical wave J 

The method according to the invention allows to determine the Reflection travel 
times of seismic events picked on 3D records of seismic data corresponding to signals 
20 picked up by seismic receivers distributed along an acquisition line (respectively along 
several acquisition lines) in response to the emission, in the subsoil, of waves from one 
or more source points, this data having been previously converted into migrated 



cylindrical wave data and interpreted (respectively in a superposition of migrated 
cylindrical wave data and interpreted). 

In cases where interpretation has been carried out on the migrated data associated 
with a single cylindrical wave, and therefore a single acquisition line, the method 
5 comprises the followin g j ^^te : 

a) defining a slowness vector {p) whose component (p x ) in a direction parallel to the 
acquisition line defines the slope of the cylindrical wave ; 

b) for a seismic receiver situated at abscissa (Xr) on the acquisition line, seeking 
0*t. abscissa (£) of the source point on the acquisition line such that a ray starting from the 

,0 s^c «*. - « on . ^ _ e^es a, i> _ poin,, wi* a 
slowness vector (p) whose component in the direction of the acquisition line is (-p x ) ; 

c) determining a travel time (t* (Xr)) by adding to the value of the travel time along the 
ray a time equal to the product of the slope of the cylindrical wave by the abscissa of the 
source point ; 

15 d) repeating stegro b) and c) for all the positions of the receivers for which a 
demigration result is wanted ; and 

e) repeating ^^^a) to d) for all the values taken by parameter (p*) for which an event 
has been picked in the corresponding migrated data. 

In cases where interpretation has been carried out on the superposition of 
20 cylindrical wave data, a superposition performed from several acquisition lines, the 
method comprises the following-! 
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a) defining a slowness vector (/?) whose component (p x ) in a direction parallel to the 
acquisition lines defines the slope of the cylindrical waves associated with the various 
acquisition lines ; 

b) for a given seismic receiver situated at abscissa (x R ) on an acquisition line, seeking 
5 0%bscissa (§) of a source point on the acquisition line such that a ray starting from the 

seismic receiver and reflecting on a picked event emerges at -said source point, with a 
slowness vector whose component in the direction of the acquisition line is (-p x ) ; 

ik 

c) determining a travel time (t e (Xr)) by adding to the value of the travel time along said 
ray a time equal to the product of the slope of the cylindrical wave by the abscissa of the 

10 source point ; 

d) repeating Sag^ b) and c) for all the positions of the receivers for which a 
demigration result is wanted ; 




e) repeating -stages a) to d) for all the acquisition lines for which a demigration result is 
wanted, 



15 f) repeating -stages a) to e) for all the values taken by parameter (px) for which an event 
has been picked in the corresponding migrated data. 

The prior processing operation leading to migrated cylindrical wave data can 
possibly involve, if necessary, AMO type preprocessing as mentioned above or a spatial 
resampling of the seismic data acquired in the field. 

20 According to an implementation mode, the travel times associated with the 

cylindrical wave events are used to implement a cylindrical-wave reflection tomography 
technique and to calculate a velocity distribution in the environment. 



According to another implementation mode, the travel times associated with the 
cylindrical wave events are converted into travel times associated with the events by 
shotpoints, by exploiting the correspondences between the cylindrical wave data and the 
shotpoint data 

According to an implementation mode, when the seismic data acquisition device is 
organized by acquisition lines, Le. when a series of source points and a series of pickups 
are arranged in one line, the synthesized shot travel times are converted to travel times 
associated with the source-receiver pairs (of the seismic data acquisition device), this 
conversion being obtained by extrapolation in the direction of the acquisition lines of 
the synthesized travel time for the source-receiver pair that is the closest to the source- 
receiver pair for which information on the travel time is wanted. 

According to an implementation mode, the seismic data being only constant- 
azimuth data, the synthesized shot travel times are converted into travel times associated 
with the source-receiver pairs, this conversion being performed by means of a double 
extrapolation, one in the direction defined by the azimuth, the second in the direction 
orthogonal to the previous direction. 

The method can be used notably to determine the velocity model by applying for 
example a reflection tomography to the travel times obtained. 

BRIEF DESCRIPTION OF THE-FfQt?RBS D^ W,r ^ 

Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of non-limitative examples, with reference 
to the accompanying drawings wherein : 
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- Figure- 1 describes the various options allowing implementation of the method 
according to the invention, and 

- FigHfe.2 compares, on a constant-offset record section, the seismic data with the 
arrival times obtained by implementing an option of the method. i 

DETAILED DESCRIPTION 0 f TH£ fdfc 

Demigration of events picked on migrated cylindrical wave data requires a specific 
technique. Wc wllkjlaccessively distinguish : 

- demigration of an event picked on migrated cylindrical wave data, 

- demigration of the superposition (stack) of migrated cylindrical wave data, 
associated with various acquisition lines, 

- obtaining kinematic data (reflection travel time) associated with source-pickup 
pairs, and 

- obtaining kinematic data associated with source-pickup pairs existing in the original 
acquisition device (i.e. for the source-pickup pairs used for migration). 

1) D emigration of an event picked on migrated cylindrical wave data 

The result of such a demigration is the time of arrival, on the records associated 
with the cylindrical wave considered (characterized by the co-ordinate along the y-axis 
of the acquisition line and by the value of parameter p* defining the slope of the 
cylindrical wave), of the reflection corresponding to the picked event. This arrival time 
tf depends on the abscissa x R of the receiver considered- Demigration is naturaUy 
conducted by seeking, for a given receiver position, the abscissa % of the source point or 
of the source points S on the line considered, such that the ray starting from the 



receiver 
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and reflecting on the picked event emerges at S with a slowness vector p whose 
component along x is -p„ the velocity model being identical to the model used for 
migration. 

Thus worded, the problem appears to be the solution of a non-linear equation 
(expressing the condition on the slowness vector) in a single unknown (abscissa § of the 
source). This equation can be solved numerically by means of any method suited for 
solution of non-linear equations, such as Newton's method for example. 

The travel time resulting from the demigration t e (x R ) will be obtained by taking the 
value of the travel time along all the rays obtained for source S, a time to which quantity 
p>£ is added. It can be noted that there is in fact an infinity of such rays, but these rays 
give, for a given arrival, the same travel time. This operation is repeated for all the 
receiver abscissas for which the demigration result is sought. The arrival times thus 
obtained constitute intrinsic information connected with the synthesized cylindrical 
wave data. They do not depend on the velocity model used for the migration (which 
justifies the fact that they can be used for determination of the velocity model), because 
the demigration operation described above is, as regards kinematics, the inverse of the 
migration of cylindrical wave data. Component p R along x of slowness vector p of the 
ray at the emergence of the seismic pickup also constitutes a quantity that is intrinsically 
linked with the data and which therefore does not depend on the velocity model used for 
the migration ; it gives the slope, on the cylindrical wave data, of the event recorded by 
the pickup. 
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2) Demigration of the superposition of migrated cylindrical wave data, associated 
with various acquisition lines 

This operation is simply carried out by repeating the operation defined in the 
previous paragraph for all the acquisition lines where a demigration result is desired and 
^delivers the same intrinsic quantities, i.e. the time of arrival of the event, on the 
cylindrical wave data associated with each line, as a function of the abscissa of the 
pickup and of the slope of this event on the x-axis, i.e. component Pr along x of the 
slowness vector when the ray emerges at the pickup. Furthermore, if we. denote the 
components on the y-axis of the slowness vector for .rays emerging at the surface 
respectively at the pickup and at the sourcej^by p y R and pi , quantity p\ + p y s also 

#A 

appears as an intrinsic quantity^JF is interpreted as the slope, on the y-axis (therefore 
from one line to the next), of the event seen for the seismic pickup abscissa considered, 
on cylindrical wave data associated with parameter p x . 

3) Obtaining kinematic data (reflection travel time) associated with source-pickup 
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lijL reflection tomography software ^jB6ying^6 process travel times associated with 
cylindrical waves is available, the data obtained in the previous paragraph is sufficient 
for updating the velocity model and for migration velocity analysis according to the 
Smart method. 

If such inversion softwares are not available, and if conventional softwares 
processing kinematic data organized by source-receiver pairs are used, the stage 
described hereafter will be carried out. 
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This stage C&isisls-m exploiting correspondences between the cylindrical wave data 
and the shotpoint data, Le. all the records obtained with the various receivers for the 
same seismic source assumed to be a point source. 

Such a correspondence can be exploited only if cylindrical wave data associated 
with several values of parameter p* defining the slope of the cylindrical waves is 
available. It (^ST be noted that migration has no role here ; the point is simply to 
transform data. The kinematic data (reflection travel time) associated with source- 
pickup pairs is obtained by proceeding acquisition line by acquisition line (the source- 
pickup pairs are assumed to be situated on these lines). 

For a given line, we have , by applying the techniques described in the previous two 
paragraphs^unction r (x R ) for various values of parameter p x . In other words, w£have 
the functioig" that we shall cati, after a notation change, r (p x , x R ), or more exactlj/a 
sampling of this function, the values taken by this functionhekig obtained for specific 
values of the (p„ Xr) pair. A function X s (p*, x R ) is defined by the relation : 

X s <J>x*x R ) = —ip x ,x R ) (1) 

a quantity that>^ cam evaluate numerically (for example by using approximations 
based on finite differences) for specific values of the (p„ Xr) pair. Formula (2) hereafter 
gives the time denoted by f (x s ,Xr) of the event at the pickup of abscissa x R in the 
shotpoint data for a shot abscissa x s = X s (p„ Xr). 



T\x s ,x R ) = t e (p x ,x R )-p x X s (p x , XR ) 
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The times thus obtained are intrinsic quantities : they are only related to the seismic 
data and they do not depend on the way migration has been implemented. In conclusion, 
it is sufficient to apply formula (2) on the various lines for all the (p^ x R ) pairs for which 

data t 6 (J)*, Xr) and ■— (p x , x R ) is available to know the arrival times of the interpreted 
VP* 

event (or, in practice, of each interpreted event) for a multitude of source-pickup pairs 
covering the surface and thus allow precise determination, via a reflection travel time 
tomography, of the velocity model. 

4) Obtaining kinematic data associated with source-pickup pairs existing in the 
original acquisition device 

Even if the data obtained by means of the technique described in the previous 
section is sufficient to update the velocity model, it may be useful to compare the times 
obtained with the non-migrated original seismic data. Such a comparison is in fact 
interesting for validating the picking (and notably the interpretation) performed on the 
migrated data and for checking the reliability of the demigration which, as mentioned 
above, involves numerical procedures based on approximations and ray tracing, a very 
sensitive tool for numerical representation of the reflectors. 

However, direct comparison is difficult since there is no reason for the source- 
pickup pairs to be those corresponding to the data used for the migration. In feet, the 
travel times provided by the method described in the previous section are obtained for 
source-pickup pairs distributed as a series of clusters distributed in 2D spaces (each one 
of these spaces corresponding to an acquisition line), whereas the original data can be 
distributed in form of a cluster distributed in a 3D space. However, even in cases where 
migration has been carried out from data organized by 3D lines, the problem remains 
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the same : there is no reason for the source-pickup pair defined by the pair of abscissa 
(Xs (px, Xr ), x R ), obtained as explained in the previous section, to appear in the original 
acquisition device. 

Interpolation or extrapolation techniques are necessary to obtain kinematic data 
5 associated with source-pickup pairs existing in the original acquisition device from the 
results obtained according to the method described in the previous section. One 
technique consists in defining a proximity criterion in a 2D or 3D space and, for a given 
source-receiver pair of the original acquisition device, in assigning thereto the travel 
time of the source-pickup pair that is the closest among those for which travel time 
10 information has been obtained by means of the method described in section 3. Such an 
approach is viable if the information extrapolation distance is not too great. In the 
opposite case, more elaborate techniques have to be implemented, such as those 
described hereafter, according to whether resampled data or constant-azimuth data is 
used. 

15 4-1 Implementation on data organi zed bv acquisition lines 



parallel to the direction defined by the azimuth, for example after spatial resampling of 
the data and/or AMO type preprocessing, as mentioned above. The source-pickup pairs 
20 involved in the implementation of the migration method are distributed as clusters 
distributed in 2D spaces, each one of these spaces corresponding to an acquisition line. 



In this case, the data is organized by acquisition line and wo want to obtain the time 



In this implementation type, wfr alao ouppos e that the seismic data is constant- 
azimuth data and that acquisition is carried out by displacing the source along lines 





of arrival of the event reflected for the source-pickup pair of the original acquisition 
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device (disp), defined by the acquisition line to which it belongs and by abscissas 
xf* and xf" . As explained above£we-haveJ extracted from all the arrival times 
calculated for the line consideree^he arrival time associated with the source-pickup pair 
defined by pair that is the closest to the time associated with pair 

among the cluster obtained in the 2D space described by pairs (x s , Xr). W 
then use formula (3) given hereafteijyto obtain the arrival time t e {xf v ,xf p ) from the 
quantities obtained during the previous processing stages : 

t e {x?> 9 xp> ) = fVt,*?) - - xfl + PA*** - xf c ) (3) 

quantity p R having been defined in sections 1, 2 above. It can be noted that this formula 
exploits the following important result : according to hypotheses 1 and 2 stated in 
section 1, for a given pickup and therefore for a given acquisition line, the slope of an 
event in the cylindrical wave data associated with parameter p x is identical to the slope 
of the event in the shotpoint collection for the shot related by relation (1) to pair (p^ Xr) 
of abscissa (x R ) of the pickup. The aforementioned extrapolation technique is of course 
not limitative. An interpolation technique can be used just as well. 

4-2 Implementation of the method on constant-azimuth data 

A 

In this implementation type, we assume that the seismic data is constant-azimuth 
data. The source-pickup pairs involved in the implementation of the migration are 
distributed as clusters distributed in a 3D space, and these pairs can be described by 
means of three parameters : the respective abscissas of the source and of the pickup, and 
their common ordinate. The problem j^nmotfl-horc in obtaining the arrival times of the 
event reflected for the source-pickup pairs of the original acquisition device (therefore 
distributed as clusters in a 3D space) from the calculated times, by using the method 
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presented in section 3. This method is then implemented for a series of fictitious 
acquisition lines of ordinate y h i representing the number of the line considered. These 
lines have to be sufficiently close to one another to allow extrapolation, on the source- 
pickup pairs of the original acquisition device, of the travel times calculated according 
to the method described in section 3. We thuo want to calculate the travel time of the 
reflected event for a source-pickup pair of the original device, this pair having 
(xp",x^ p ) as the abscissa and y 6 * as the common ordinate. 

We t h erefor^^ estEaetr from all the arrival times calculated for the series of lines 
mentioned above, the time associated with the source^pickup pair defined by pair 
xf c ,x^) and by ordinate y, that is the closes^ according to the metric defined in the 



3D space, to the time associated with the source-pickup pair considered of the original 
device. We-then-«se formula (4) hereafte^jto obtain the arrival time t e (xf sp ,xf sp ^y disp ) 
from the quantities obtained during the previous processing stages : 

axp> 9 xp> 9 y**>) = f{x<t,xf\ yi ) - pA4* - *f c ) + PxOfr -xf°) + iP y R + pSW** - y t ) 

Here again, the extrapolation technique is given by way of non limitative example. 
15 An interpolation technique can be used just as well. 

Figore 2 compares, on a constant-offset section (1500 m here) obtained after 3D 
seismic prospecting operations in the North Sea, the seismic data with the arrival times 
obtained by implementing the technique described above. The reflection can be clearly 
seen at the base of the salt structure. 
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CLAIMS 

1) A method for determining reflection travel times of seismic events picked on 3D 
records of seismic data corresponding to signals picked up by seismic receivers 
distributed along an acquisition line in response to the emission in the subsoil of waves 
5 from source points, 313s data having been previously converted to cylindrical wave data 
and interpreted, c ^ackfil^iin^ it ajiiiyiises the fulluwhig mage s : 

a) defining a slowness vector (p) whose component (p x ) in a direction parallel to the 
acquisition line defines the slope of the cylindrical wave ; 

b) for a seismic receiver situated at^abscissa (x R ) on the acquisition line, seeking 
10 ^abscissa (£) of the source point on the acquisition line such that a ray starting from the 

seismic receiver and reflecting on a picked event emerges at the source point, with a 
slowness vector (p) whose component in the direction of the acquisition line is (-p x ) ; 

c) determining a travel time (t* (x R )) by adding to the value of the travel time along jJkI 
ray a time equal to the product of the slope of the cylindrical wave by the abscissa of the 

15 source point; 




tor all the positions of th 

demigration result is wanted ; and 



d) repeating -stages b) and c) for all the positions of the ^receivers for which a 




e) repeating «agef5 a) to d) for all the values taken by parameter (p x ) for which an event 
has been picked in the corresponding migrated data. 

20 2) A method for determining reflection travel times of seismic events picked on 3D 

records of seismic data corresponding to signals picked up by seismic receivers 
distributed along several acquisition lines in response to the emission in the subsoil of 
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waves from one or more source points, this data having been previously converted to a 



superposition of migrated cylindrical wave data and interpreted, cnaraotorizo c bL lliflr it 
ge mpil^es ( lie lolk>Wiiig"Siage s : 

a) defining a slowness vector (p) whose component (p x ) in a direction parallel to the 
5 acquisition lines defines the slope of the cylindrical waves associated with the various 

acquisition lines ; 

b) for a given seismic receiver situated at an abscissa (x R ) on an acquisition line, seeking 
Ql* abscissa <£) of a source point on the acquisition line such that a ray starting from the 

seismic receiver and reflecting on a picked event emerges at^afd source point, with a 
10 slowness vector (p) whose component in the direction of the acquisition line is (-p x ) ; 

tL 

c) determining a travel time (t* (x R )) by adding to the value of the travel time along said 
ray a time equal to the product of the slope of the cylindrical wave by the abscissa of the 
source point ; 




d) repeating stages b) and c) for all the positions of the receivers for which a 
15 demigration result is wanted ; 

e) repeating stages a) to d) for all the acquisition lines for which a demigration result is 



eg) A 



want 

f) repeating staged a) to e) for all the values taken by parameter (p x ) for which an event 
has been picked in the corresponding migrated data. 

20 3) A method as claimed in^giy one o^clain^ 1 ch^actoizod in that the travel 

times associated with the cylindrical wave events are used to implement a cylindrical- 
wave reflection tomography technique and to calculate a velocity distribution in the 
environment. 
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4) A method as claimed in any &m of clainy I-w2, c baractgnzed m th-rf the travel 
times associated with the cylindrical wave events are converted to travel times 
associated with shotpoint events, by exploiting correspondences between the cylindrical 
wave data and the shotpoint data. 

5 5) A method as claimed in claim 4, characterized in that , the seismic data b^mg 

organized by acquisition lines, the synthesized shot travel times are converted to travel 
times associated with source-receiver pairs, tltft conversion being carried out by 
O extrapolation, in the direction of the acquisition lines, of the travel time synthesized for 

t£6 source-receiver pair that is tjfe closest to the source-receiver pair of the acquisition 
10 device for which information on the travel time is wanted. 

6) A method as claimed in claim 4, c haracterized in that , the seismic data bpfiig 



r : 



constant-azimuth data, the synthesized shot travel times are converted to travel times 
C associated with the source-receiver pairs, conversion being carried out by means of 

a double extrapolation, the first one in the direction defined by the azimuth, the second 
15 in the direction orthogonal to the previous direction 

7) Application of the method as claimed in any one of clainyS 4 tor 6 fm determining 
velocity model by applying, -for exampl e, a reflection tomography to ip& travel times 
obtained. 
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MSTITUT FRANgAIS DU PETROLE 

METHOD INTENDED TO OBTAIN REFLECTION TRAVEL TIMES FROM 
AN INTERPRETATION OF MIGRATED CYLINDRICAL WAVE 

SEISMIC DATA 

Inventors : Patrick LAILLY and Frederic ASSOULINE 
ABSTRACT 

f Method intended to obtain reflection travel times from an interpretation of seismic 
data in migrated cylindrical waves, for a given value of the parameter defining the slope 
of these waves, or the superposition of such data associated with various substantially 
parallel acquisition lines, this parameter possibly taking successively several values. 

/^The method comprises the fu^|vi^g slagus : a) defining a slowness vector {p) 
^hose component (p x ) in a direction parallel to the acquisition lines defines the slope of 
the cylindrical wavj^; b) for a given position of a seismic receiver of abscissa (x R ) on an 
acquisition line, seeldng^scissa <£) of the source^such that a ray starting from the 
seismic receiver and reflecting on a picked event emerges at the source on the 
acquisition line considered, with a slowness vector (J?) whose component in the 
direction parallel to the acquisition lines is (-pj; c) determining a travel time (t e (x^jby 
adding to the value of the travel time along said ray a time equal to the product of slope 
(p x ) of the cylindrical wave by the abscissa of the source poi^; d) repeating iia^s b) 
and c) for all the positions of the receivers for which a demigration result is wanted ; 
and e) repeating Qte^ a) to d) for all the acquisition lines for which a demigration 
result is wanted and for all the values taken by parameter (p x ). 

Applications : migration velocity analysis for determination of the propagation 
velocities distribution of seismic waves in the subsoil for example^ 
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METHOD INTENDED TO OBTAIN REFLECTION TRAVEL TIMES FROM 
AN INTERPRETATION OF MIGRATED CYLINDRICAL WAVE 

SEISMIC DATA 



Inventors : Patrick LAILLY and Frederic ASSOULINE 



BACKGROUND OF THE INVENTION 

Field Of The Invention 

The present invention relates to a method for obtaining reflection travel times from 
an interpretation of seismic data in migrated cylindrical waves, for a given value of a 
parameter defining the slope of these waves, or from the superposition of such data 
associated with various substantially parallel acquisition lines. 



Description of the Prior Art 

Demigration of seismic data is the inverse of migration, which is a conventional 
seismic data imaging method essentially consisting, knowing the value of a wavefield at a 
known depth, for example at the surface, as well as a distribution model of the 
propagation velocities of the waves in the subsoil, in modeling the propagation of the 
source field and the backpropagation of the reflection data recorded, and in seeking phase 
coherences between these two modeled fields. Migration is particularly useful for 
interpretation of seismic data acquired on complex structures. 
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It also allows to access, via demigration, the arrival times of the reflections 
associated with the picked events. This operation is notably carried out for determining, 
by means of kinematic methods such as migration velocity analysis methods, the velocity 
distribution in the subsurface. This determination constitutes a key stage in complex 
structures imaging. An example of such a method is described by: 

Stork, C, 1992; «Reflection Tomography in the Postmigrated Domain»; Geophysics, 
57, 680-692. 

A particularly effective migration velocity analysis method, referred to as Smart 
(«Sequential Migration Aided Reflection Tomography ») is notably described by: 
Ehinger, A. and Lailly, P., 7995, Velocity Model Determination by the SMART 
Method, Part 1: Theory: 65th Ann. Internal Mtg., Soc. Expl. Geophys., Expanded 
Abstracts, pp. 739-742. 

Reflection tomography is used, in such methods, for determining the velocity 
model The reflection travel times (reflection tomography data) are obtained, in the 
complex parts, by demigrating the picked events upon interpretation of the migrated 
data. This method has been successfully implemented for 3D imaging of complex 
structures by using 2D prestack migration algorithms (and the associated 
demigration) on a set of seismic profiles extracted from a 3D acquisition. However, 
such an approach is efficient only in a very specific context: 

- on the geologic plane, the structure must vary relatively slowly in a direction 
(referred to as « strike» in English literature); 

- the seismic profiles must be acquired in the orthogonal direction (therefore the 



complexity direction) referred to as «dip» in the English literature. 

On the other hand, this approach is difficult to use with complex structures 
without a privileged direction: such structures require prestack 3D migration 
techniques, which involve long calculation times, especially if imaging is carried out 
on a volume. 

French Patent 2,784,195, filed by the assignee, notably describes another 
seismic data migration method known in the art as Alchemig. This method allows 
performing a 3D prestack migration of seismic events and imaging of the volumes of 
an underground zone, with an attractive calculation time, from a series of a number 
Ns of reflection survey cycles. 

It comprises successive emission of elementary wavefields, each defined by the 
association of a seismic signal W(t) and of a determined emission point from a series 
of emission points S with 1 < i < Ns, reception, by seismic receivers placed at 
positions Rf , of the seismic signals reflected by the zone in response to each one of 
these wavefields, and recording of the various signals received by each seismic 
receiver in form of time-dependent seismic traces d/ (t). For a given velocity model, 

a) a slowness vector p (homogeneous to the inverse of a velocity) whose two 
components p y and p x can take a previously defined series of values is defined, 

b) a time lag function to( p J) is defined for a given slowness vector p and for a 
given emission point S i, 



c) a time lag function to( p A) is applied to each elementary wavefield associated 

with emission point S i, and a composite wavefield is formed at the surface by 
spatiotemporal superposition of the various elementary wavefields to which such 
a time lag is applied, 

d) a time lag to(pJ) is applied to each seismic trace d\ (t) marked by pair (i,j) and a 

composite trace field is formed at the surface by spatiotemporal superposition of 
the various seismic traces to which such a time lag is applied, 

e) a migration of the composite trace field is carried out by using as the wavefield 
the composite wavefield, which is done by modeling the propagation of the 
composite wavefield and the backpropagation of the composite trace field, and 
by suitably combining the two composite fields thus modeled at any point of the 
zone to be imaged, 

f) steps (c to e) are repeated for all the values taken by components p y and p x of 
vector p , and 

g) for any fixed value of the second component p x of the vector p , the results of 
these various combinations are summed so as to obtain a migrated image 
associated with this fixed value of p x , thus performing a prestack migration. 

It is assumed that, even if the technique proves robust in case of violation of these 
hypotheses, the seismic data is constant-azimuth data (the azimuth being defined as the 
direction of the source-pickup bipoints) (hypothesis 1) and that acquisition is carried out 
by displacing the source along lines parallel to the direction defined by the azimuth 
(hypothesis 2). 



If this is not the case, it can be remedied by first carrying out a preprocessing 
operation, for example according to the technique referred to as AMO (« Azimuth 
Moveout »), well-known to the man skilled in the art, to adjust the data obtained from 
standard marine acquisitions. 

The co-ordinates system selected here is such that the x-axis is parallel to the 
direction of the acquisition lines and the y-axis thus represents the position of an 
acquisition line. Hereafter the measurement by cylindrical wave data is defined, along an 
acquisition line, of the seismic response of the subsurface to a cylindrical wave (wave 
generated by a source line and whose phase varies linearly with the abscissa along the 
acquisition line) whose axis coincides with the acquisition line considered. Considering 
this definition, the Alchemig technique allows obtaining the superposition, along the 
various acquisition lines, of migrated cylindrical wave data, these cylindrical waves being 
associated with an acquisition line and with a predetermined value of parameter p x 
defining the slope of the cylindrical wave. The procedure consists in applying to the data 
linear phase shifts parametered by vector p whose component along the profiles axis is 
p x , and in calculating such a superposition by applying to the phase shifted data a plane- 
wave migration software. Postmigration stacking is obtained by summing the results 
obtained for the different values taken by p x . We suppose that, by suitable selection of the 
origin of the mark or of the phase shift function, the point of the surface where the phase 
shifts are zero has (x =0, y =0) as the co-ordinates. 



It can be noted that the technique described above is one way among many to obtain 
superposition of migrated cylindrical wave data associated with parallel acquisition lines 
and with a given value of parameter p x . Another procedure, can for example, is described 
by Claerbout, J.R, 1971, Towards a Unified Theory of Reflector Mapping; in 
Geophysics, 36, No. 3, 467-481, in calculating the result by carrying out, for each 
acquisition line: 

- propagation of the cylindrical wave and back propagation of the cylindrical 
wave data, and 

- calculation, at any point of the space, of the cross correlation of the source 
field whose propagation has been modeled and of the back propagated data 
field, which allows obtaining, for the acquisition line considered, the migrated 
cylindrical wave data; 

and by stacking (superposing) the migrated cylindrical wave data obtained for 
the various acquisition lines. 

If, for several values of parameter p X9 the superposition of migrated cylindrical wave 
data is associated with various acquisition lines, it is possible, as explained in the 
aforementioned patent, to control the quality of the velocity model used for the migration 
and also to update the velocity model by exploiting the deformations of the events when 
going from one cylindrical wave parameter to the next (migration velocity analysis). 

For updating, according to the aforementioned Smart method (or others), the events 
have to be picked, for various values of parameter p X9 in the superposition of the migrated 



cylindrical wave data, a superposition corresponding to the various acquisition lines, then 
the picked events have to be demigrated and the model is updated by reflection 
tomography. Even it is required to be known how to demigrate the interpreted migrated 
events when interpretation has been performed on migrated cylindrical wave data or the 
superposition obtained for various acquisition lines. 

SUMMARY OF THE INVENTION 

The method according to the invention allows determination of the reflection travel 
times of seismic events picked on 3D records of seismic data corresponding to signals 
picked up by seismic receivers distributed along an acquisition line (respectively along 
several acquisition lines) in response to the emission, in the subsoil, of waves from one or 
more source points, this data having been previously converted into migrated cylindrical 
wave data and interpreted (respectively in a superposition of migrated cylindrical wave 
data and interpreted). 

In cases where interpretation has been carried out on the migrated data associated 
with a single cylindrical wave, and therefore a single acquisition line, the method 
comprises the following steps: 

a) defining a slowness vector ( p ) whose component (p x ) in a direction parallel to 

the acquisition line defines the slope of the cylindrical wave; 
b) for a seismic receiver situated at abscissa (Xr) on the acquisition line, seeking an 
abscissa (f) of the source point on the acquisition line such that a ray starting 
from the seismic receiver and reflecting on a picked event emerges at the source 



point, with a slowness vector ( p ) whose component in the direction of the 
acquisition line is (-p x )l 

c) determining a travel time (t e (X R )) by adding to the value of the travel time along 
the ray a time equal to the product of the slope of the cylindrical wave by the 
abscissa of the source point; 

d) repeating steps (b and c) for all the positions of the receivers for which a 
demigration result is wanted ; and 

e) repeating steps (a to d) for all the values taken by parameter (p x ) for which an 
event has been picked in the corresponding migrated data. 

In cases where interpretation has been carried out on the superposition of cylindrical 
wave data, a superposition performed from several acquisition lines, the method 
comprises the following steps: 

a) defining a slowness vector (p ) whose component (p x ) in a direction parallel 
to the acquisition lines defines the slope of the cylindrical waves associated 
with the various acquisition lines; 

b) for a given seismic receiver situated at abscissa (X R ) on an acquisition line, 
seeking an abscissa (f) of a source point on the acquisition line such that a ray 
starting from the seismic receiver and reflecting on a picked event emerges at 
the source point, with a slowness vector ( p ) whose component in the direction 
of the acquisition line is (-p x ); 

c) determining a travel time (t e (Xr)) by adding to the value of the travel time 
along the ray a time equal to the product of the slope of the cylindrical wave by 
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the abscissa of the source point; 

d) repeating steps (b and c) for all the positions of the receivers for which a 
demigration result is wanted; 

e) repeating steps (a to d) for all the acquisition lines for which a demigration 
result is wanted; 

f) repeating steps (a to e) for all the values taken by parameter (p x ) for which an 
event has been picked in the corresponding migrated data. 

The prior processing operation leading to migrated cylindrical wave data can possibly 
involve, if necessary, AMO type preprocessing as mentioned above or a spatial 
resampling of the seismic data acquired in the field. 

According to an implementation mode, the travel times associated with the 
cylindrical wave events are used to implement a cylindrical-wave reflection tomography 
technique and to calculate a velocity distribution in the environment. 

According to another implementation mode, the travel times associated with the 
cylindrical wave events are converted into travel times associated with the events by 
shotpoints, by exploiting the correspondences between the cylindrical wave data and the 
shotpoint data. 

According to an implementation mode, when the seismic data acquisition device is 
organized by acquisition lines, i.e. when a series of source points and a series of pickups 



are arranged in one line, the synthesized shot travel times are converted to travel times 
associated with the source-receiver pairs (of the seismic data acquisition device), this 
conversion being obtained by extrapolation in the direction of the acquisition lines of the 
synthesized travel time for the source-receiver pair that is the closest to the source- 
receiver pair for which information on the travel time is wanted. 

According to an implementation mode, the seismic data being only constant-azimuth 
data, the synthesized shot travel times are converted into travel times associated with the 
source-receiver pairs, this conversion being performed by means of a double 
extrapolation, one in the direction defined by the azimuth, the second in the direction 
orthogonal to the previous direction. 

The method can be used notably to determine the velocity model by applying for 
example a reflection tomography to the travel times obtained. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of non-limitative examples, with reference to 
the accompanying drawings wherein: 

Fig.l describes the various options allowing implementation of the method according to 
the invention, and 

Fig.2 compares, on a constant-offset record section, the seismic data with the arrival 
times obtained by implementing an option of the method. 



10 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Demigration of events picked on migrated cylindrical wave data requires a 

specific technique. Successively distinguished are: 

- demigration of an event picked on migrated cylindrical wave data, 
-demigration of the superposition (stack) of migrated cylindrical wave data, 
associated with various acquisition lines, 

- obtaining kinematic data (reflection travel time) associated with source-pickup 
pairs, and 

- obtaining kinematic data associated with source-pickup pairs existing in the 
original acquisition device (i.e. for the source-pickup pairs used for migration). 

1) Dem igration of an event picked on migrated cylindrical wave data 

The result of such a demigration is the time of arrival, on the records associated with 
the cylindrical wave considered (characterized by the co-ordinate along the y-axis of the 
acquisition line and by the value of parameter p x defining the slope of the cylindrical 
wave), of the reflection corresponding to the picked event. This arrival time t e depends on 
the abscissa X R of the receiver considered. Demigration is naturally conducted by 
seeking, for a given receiver position, the abscissa £ of the source point or of the source 
points S on the line considered, such that the ray starting from the receiver and reflecting 
on the picked event emerges at S with a slowness vector p whose component along x is 
-Px, the velocity model being identical to the model used for migration. 
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Thus worded, the problem appears to be the solution of a non-linear equation 
(expressing the condition on the slowness vector) in a single unknown (abscissa £ of the 
source). This equation can be solved numerically by means of any method suited for 
solution of non-linear equations, such as Newton's method for example. 

The travel time resulting from the demigration t e (X R ) will be obtained by taking the 
value of the travel time along all the rays obtained for source S, a time to which quantity 
p x % is added. It can be noted that there is in fact an infinity of such rays, but these rays 
give, for a given arrival, the same travel time. This operation is repeated for all the 
receiver abscissas for which the demigration result is sought. The arrival times thus 
obtained constitute intrinsic information connected with the synthesized cylindrical wave 
data. They do not depend on the velocity model used for the migration (which justifies the 
fact that they can be used for determination of the velocity model), because the 
demigration operation described above is, as regards kinematics, the inverse of the 
migration of cylindrical wave data. Component p R along x of slowness vector p of the 
ray at the emergence of the seismic pickup also constitutes a quantity that is intrinsically 
linked with the data and which therefore does not depend on the velocity model used for 
the migration; it gives the slope, on the cylindrical wave data, of the event recorded by the 
pickup. 

2) Demigration of the superposition of migrated cylindrical wave data, associated 

with various acquisition lines 
This operation is simply carried out by repeating the operation defined in the 
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previous paragraph for all the acquisition lines where a demigration result is desired and 
delivers the same intrinsic quantities, i.e. the time of arrival of the event, on the 
cylindrical wave data associated with each line, as a function of the abscissa of the pickup 
and of the slope of this event on the x-axis, i.e. component /?r along x of the slowness 
vector when the ray emerges at the pickup. Furthermore, if the components on the y-axis 
of the slowness vector for rays emerging at the surface respectively at the pickup and at 
the source are denoted by p y R and p y , quantity p y R + p y s also appears as an intrinsic 
quantity and is interpreted as the slope, on the y-axis (therefore from one line to the next), 
of the event seen for the seismic pickup abscissa considered, on cylindrical wave data 
associated with parameter p x . 

3) Obtaining kinematic data (reflectiontravel time) associated with source-pickup pairs 
If reflection tomography software which processes travel times associated with 
cylindrical waves is available, the data obtained in the previous paragraph is sufficient 
for updating the velocity model and for migration velocity analysis according to the Smart 
method. 

If such inversion softwares are not available, and if conventional softwares processing 
kinematic data organized by source-receiver pairs are used, the steps described hereafter 
will be carried out. 

This step exploits correspondences between the cylindrical wave data and the 
shotpoint data, i.e. all the records obtained with the various receivers for the same seismic 
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source assumed to be a point source. 



Such a correspondence can be exploited only if cylindrical wave data associated with 
several values of parameter p x defining the slope of the cylindrical waves is available. It 
should be noted that migration has no role here: the point is simply to transform data. The 
kinematic data (reflection travel time) associated with source-pickup pairs is obtained by 
proceeding acquisition line by acquisition line (the source-pickup pairs are assumed to be 
situated on these lines). 

For a given line, it is obtained by applying the techniques described in the previous 
two paragraphs, a function t e (X R ) for various values of parameter p x . In other words, the 
function is obtained that is named after a notation change, t e (p x , X R ), or more exactly for 
a sampling of this function, the values taken by this function are obtained for specific 
values of the (p x , X R ) pair. A function X s (p x , X R ) is defined by the relation: 

dt e 

X s(P x >*r) = —(P x ,x r ) (1) 

a quantity that can be evaluated numerically (for example by using approximations based 
on finite differences) for specific values of the (p x , X R) pair. Formula (2) hereafter gives 
the time denoted by t e (Xs,X R ) of the event at the pickup of abscissa X R in the shotpoint 
data for a shot abscissa X s = X s (p x , X R ). 

T m (x s ,x R ) = t e (j> x9 x R )- Px X s (j> x9 x R ) (2) 

The times thus obtained are intrinsic quantities: they are only related to the seismic 
data and they do not depend on the way migration has been implemented. In conclusion, 
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it is sufficient to apply formula (2) on the various lines for all the (p x , x R ) pairs for which 
data t e (p x , x R ) and "^~ (P*> x *) is available to know the arrival times of the interpreted 
event (or, in practice, of each interpreted event) for a multitude of source-pickup pairs 
covering the surface and thus allow precise determination, via a reflection travel time 
tomography, of the velocity model. 

4) Obtaining kinematic data a ssociated with source-pickup pairs existing in the original 
acquisition device 

Even if the data obtained by means of the technique described in the previous section is 
sufficient to update the velocity model, it may be useful to compare the times obtained 
with the non-migrated original seismic data. Such a comparison is in fact interesting for 
validating the picking (and notably the interpretation) performed on the migrated data and 
for checking the reliability of the demigration which, as mentioned above, involves 
numerical procedures based on approximations and ray tracing, a very sensitive tool for 
numerical representation of the reflectors. 

However, direct comparison is difficult since there is no reason for the source-pickup 
pairs to be those corresponding to the data used for the migration. In fact, the travel times 
provided by the method described in the previous section are obtained for source-pickup 
pairs distributed as a series of clusters distributed in 2D spaces (each one of these spaces 
corresponding to an acquisition line), whereas the original data can be distributed in form 
of a cluster distributed in a 3D space. However, even in cases where migration has been 
carried out from data organized by 3D lines, the problem remains the same: there is no 
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reason for the source-pickup pair defined by the pair of abscissa (X s (p x , X R ), X R ), 
obtained as explained in the previous section, to appear in the original acquisition device. 

Interpolation or extrapolation techniques are necessary to obtain kinematic data 
associated with source-pickup pairs existing in the original acquisition device from the 
results obtained according to the method described in the previous section. One technique 
consists in defining a proximity criterion in a 2D or 3D space and, for a given source- 
receiver pair of the original acquisition device, in assigning thereto the travel time of the 
source-pickup pair that is the closest among those for which travel time information has 
been obtained by means of the method described in section 3. Such an approach is viable 
if the information extrapolation distance is not too great. In the opposite case, more 
elaborate techniques have to be implemented, such as those described hereafter, 
according to whether resampled data or constant-azimuth data is used. 

4-1 Implementation on data organized by acquisition lines 

In this implementation type, it is assumed that the seismic data is constant-azimuth 
data and that acquisition is carried out by displacing the source along lines parallel to the 
direction defined by the azimuth, for example after spatial resampling of the data and/or 
AMO type preprocessing, as mentioned above. The source-pickup pairs involved in the 
implementation of the migration method are distributed as clusters distributed in 21) 
spaces, each one of these spaces corresponding to an acquisition line. 

In this case, the data is organized by acquisition line and it is desired to obtain the 
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time of arrival of the event reflected for the source-pickup pair of the original acquisition 
device (disp), defined by the acquisition line to which it belongs and by abscissas xf p 
and xf p . As explained above, extracted from all the arrival times calculated for the line 
considered is the arrival time associated with the source-pickup pair defined by pair (xf c , 
xf c ) that is the closest to the time associated with pair (xf p , xf p ) among the cluster 
obtained in the 2D space described by pairs (X s , Xr). The formula (3) given hereafter is 
used to obtain the arrival time f(xf p , Xf p ) from the quantities obtained during the 
previous processing stages: 

fc*F,4*) - - pAx p> _ + Mx p> _ xt) (3) 

quantity ^ R having been defined in sections 1, 2 above. It can be noted that this formula 
exploits the following important result: according to hypotheses 1 and 2 stated in section 
1, for a given pickup and therefore for a given acquisition line, the slope of an event in the 
cylindrical wave data associated with parameter p x is identical to the slope of the event in 
the shotpoint collection for the shot related by relation (1) to pair (p x , x R ) of abscissa 
(x R ) of the pickup. The aforementioned extrapolation technique is of course not limitative. 
An interpolation technique can be used just as well. 

4-2 Implementation of the method on constant-azimuth data 

In this implementation type, it is assumed that the seismic data is constant-azimuth 
data. The source-pickup pairs involved in the implementation of the migration are 
distributed as clusters distributed in a 3D space, and these pairs can be described by 
means of three parameters the respective abscissas of the source and of the pickup, and 
their common ordinate. The problem is in obtaining the arrival times of the event 
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reflected for the source-pickup pairs of the original acquisition device (therefore 
distributed as clusters in a 3D space) from the calculated times, by using the method 
presented in section 3. This method is then implemented for a series of fictitious 
acquisition lines of ordinate y i? i representing the number of the line considered. These 
lines have to be sufficiently close to one another to allow extrapolation, on the source- 
pickup pairs of the original acquisition travel times calculated according to the method 
described in section 3. It is desired to calculate the travel time of the reflected event for a 
source-pickup pair of the original device, this pair having (xf p , xf p ) as the abscissa and 
y dlsp as the common ordinate. 

From all the arrival times calculated for the series of lines mentioned above, the time 
associated with the source-pickup pair defined by pair {xf c t x™ lc ) and by ordinate y, that 
is the closest is extracted according to the metric defined in the 3D space, to the time 
associated with the source-pickup pair considered of the original device. The formula (4) 
hereafter is used to obtain the arrival time f(x f p ,xf p , y disp ) from the quantities obtained 
during the previous processing stages: 

axr.xl?^) = t\xr,xt, yi )- pAxf - xf) + p K Ofr - xf) + + /*X/* 

Here again, the extrapolation technique is given by way of non limitative example. An 
interpolation technique can be used just as well. 

Fig. 2 compares, on a constant-offset section (1500 m here) obtained after 3D seismic 
prospecting operations in the North Sea, the seismic data with the arrival times obtained 
by implementing the technique described above. The reflection can be clearly seen at the 
base of the salt structure. 
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ABSTRACT 

A method intended to obtain reflection travel times from an interpretation of seismic data 
in migrated cylindrical waves, for a given value of the parameter defining the slope of these 
waves, or the superposition of such data associated with various substantially parallel 
acquisition lines, this parameter possibly taking successively several values. The method 
comprises the steps: 

a) defining a slowness vector ( p ) 

b) for a given position of a seismic receiver of abscissa (X R ) on an 
acquisition line, seeking, the abscissa (f) of the source; 

c) determining a travel time (t e (Xr) ) 

d) repeating steps (b and c) for all the positions of the receivers for which 
demigration result is wanted; and 

e) repeating steps (a to d) for all the acquisition lines for which a demigration result is 
wanted and for all the values taken by parameter (p x ). 
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